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ABSTRACT 



Aims. For the first time, the motion of granules (solar plasma on the surface on scales larger than 2.5 Mm) has been followed over the 

entire visible surface of the Sun, using SDO/HMI white-light data. 

Methods. Horizontal velocity fields are derived from image correlation tracking using a new version of the coherent structure tracking 

algorithm.The spatial and temporal resolutions of the horizontal velocity map are 2.5 Mm and 30 min respectively . 

Results. From this reconstruction, using the multi-resolution analysis, one can obtain to the velocity field at diff'erent scales with its 

derivatives such as the horizontal divergence or the vertical component of the vorticity. The intrinsic error on the velocity is ~ 0.25 

km s~^ for a time sequence of 30 minutes and a mesh size of 2.5 Mm.This is acceptable compared to the granule velocities, which 

range between 0.3 km s~^ and 1.8 km s~^ A high correlation between velocities computed from Hinode and SDO/HMI has been found 

(85%). From the data we derive the power spectrum of the supergranulation horizontal velocity field, the solar diff'erential rotation, 

and the meridional velocity. 

Key words. The Sun: Atmosphere - The Sun: Granulation 



1. Introduction 

The Sun is a star whose plasma flows are the main source of its 
magnetic evolution and activity. The coupling of these motions 
and the magnetic field in the convection zone is the driver of the 

] magnetic activity through the dynamo process. The description 
of the physical properties of the convective zone requires the 
knowledge of these motions at all scales in space and time. 

Determining plasma motions inside the Sun is a difl&cult task. 
While local helioseismology can be used to infer vector flows in 

] three dimensions just below the solar surface, this is only pos- 
sible for spatial scales larg er than 5 Mm an d tem poral scales 
longe r than several hour s (Ijackiewicz et al. i 2008; Gizon et aP 
l2QlQl:ISvanda et al.ll20Tlh . 

From the analysis of the proper motion of photos pheric 
struct ures (solar granules), it has been shown (Rieutord et al.l 
l200lh that it is possible to get the horizontal plasma flow on 
the Sun surface. In mo re detail, using numer ical simulation in 
a large horizontal box, iRieutord et al.l (12001 1) show that gran- 
ule motions are highly correlated with horizontal flows when the 
scale is larger than ~ 2500 km; at smaller scale, granule motions 
should be considered as (solar) turbulent noise. Such techniques 
have been used on relatively small fields of view (few arcmin- 
utes) and usually located at the disk center. 
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Very recently, Hinode and SDO satellites produced time se- 
quences of images of the solar surface without disturbance from 
the terrestrial atmosphere (seeing). These observations allow us 
to simultaneously measure horizontal velocities at diflTerent spa- 
tial resolutions and to study long-time sequences of velocities 
over the fufl Sun (HMI/SDO). 

In this paper, we describe a method for determining hori- 
zontal velocities from proper motions of the solar structures ob- 
served on the full-disk Sun. For the first time, the motion of 
granules (solar plasma flows on the surface on scales larger than 
2.5 Mm) has been followed over the full visible Sun surface. 
Horizontal velocity fields are derived from granule tracking us- 
ing a new version of the co herent structure tracking (GST) algo- 
rithm (IRieutord et al.ll2007b . 

We first present in some detail the algorithm based on gran- 
ule tracking, which is able to give a reconstruction of the velocity 
field at all scales larger than the sampling scale. The GST also 
ofifers the possibility of selecting specific structures according 
to their nature, size, lifetime, etc. and study their motion. In the 
next section we discuss the diflTerent steps of the algorithm, the 
segmentation and interpolation processes, and a comparison of 
Hinode and SDO flow maps. We illustrate the method with appli- 
cations such as the calculation of the kinetic power spectrum of 
supergranulation, and the measurement of the solar difiTerential 
rotation and meridional flow at the central meridian. Discussion 
and conclusions follow. 
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Fig. 1. Full Sun HMI/SDO white-light on August 30, 2010 (left) and the segmented map, in which about 500 000 granules are 
detected (right). 



2. Observations 

2.1. Hinode observations 

We used data sets of the Solar Opt ical Telescope ( SOT) , 
on-board the Hino d^ mission (e.g. ISuematsu et al.l 120081: 
llchimoto et al.ll2004l) . The SOT has a 50 cm primary mirror with 
a spatial resolution of about 0.2 ''at 550 nm. For our study, 
we used blue continuum observations at 450.45 nm from the 
Hinode/SOT BFI (BroadBand Filter Imager). The observations 
were recorded on August 30 2010 from 08:04:33 to 10:59:36 
UT. Solar rotation is compensated to observe exactly the same 
region of the Sun. The time step is 60 sec and the field of view 
is 76.7 '' X 76.7 '' with a pixel scale of (Y! 1089. After alignment, 
the useful field-of- view is reduced to 76V0 x 74V9. 

To remove the eff'ects of the oscillations, we applied a sub- 
sonic Fourier filter. This filter is defined by a cone in the k - co 
space, where k and oj are spatial and temporal frequencies. All 
Fourier components with oj/k < V = 7kms~^ were retained to 
keep only convective motions (JTitle et al.lll989h . 

2.2. SDO-HMI observations 

The Helioseismic and Magnetic Imager onboard the Solar 
Dynamics Observatory (HMI/SDO) yields uninterrupted high- 
resolution observations of the entire disk. This gives a unique 
opportunity for mapping surface flows on various scales (spa- 
tial and temporal). Using the HMI/SDO white-light data from 
August 30 2010, we derived horizontal velocity fields from im- 
age granulation tracking using a newly developed version of the 
CST algorithm. The time step is 45 seconds with a pixel scale 
of a'5. HMI/SDO white-light data from April 10 2010, July 10 



^ The Hinode spacecraft launched in 2006, was designed and is 
now operated by JAXA (Japanese Space Authorities) in cooperation 
with NASA (National Aeronautics and Space Administration) and ESA 
(European Space Agency). 



2010 and May 12,13 201 1 are used to measure solar differential 
rotation during 24h for each of them. 



3. The new CST method 

Granule tracking is a difficult task because of the complex evolu- 
tion on granules depending of their size. Granule fragmentation 
and mixing must be managed very carefully to avoid the gener- 
ating noise. We pursue here the development of a method called 
coherent structure tracking, or CST, which deter mines the hor- 
izontal motion of granules in the field of view (Rieutord et al.l 
2007). The preceding version was developed for a field of a few 
arcminutes, generally located at the disk center of the Sun where 
the solar rotation was removed by alignment of the frames. The 
application of the CST to HMI/SDO data requires a new granule 
time labeling methode which takes care of the motion caused 
by the solar rotation to avoid misidentification. Indeed, in the 
preceding CST version the time labeling was processed by fol- 
lowing the bary center of the granule. The solar rotation that is 
present in HMI/SDO data and the evolution of granules like the 
birth of a new (small) granule close to an existing granule be- 
tween two frames can lead to misidentification of bary centers. 
This leads to a bad temporal labeling and generates noticeable 
noise in the final velocity maps. 



3.1. Segmentation and granule identification 

To identify a granule one needs to establish a local criterion to 
decide whether a given pixel belongs to a given structure. This 
criterion needs to be local to avoid threshold effects due to large- 
scale variations of the intensity. The most efficient segmentation 
algorithm for solar granulation is an adaptated versio n of St rous' 
segmentati on (curvature-based c riterion Strous L.j (Il994l) ) de- 
scribed in iRieutord et al.l (l2007l) . This image segmentation is 
very efficient at removing large-scale intensity fluctuations. An 
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Fig. 2. Close-up view of the field around a spot as delineated in Fig. [T] (top left) with its segmentation (top right). Below, resulting 
velocity field from granule tracking (note the clear signature of the moat flow associated with the spot near x=50,y=100). 



example is shown in Fig.[T] To summarize, it consists of the fol- 
lowing steps: 

- Calculation of the "minimal curvature image": for each 
pixel, the minimal curvature among the four directions is 
computed. 

- Detection of the granules as non-negative curvature pixels in 
the minimal curvature image. 

- Extension of the detected granules with points whose mini- 
mal curvature value is higher than a given (negative) thresh- 
old, while keeping a minimal distance of one pixel between 
each granule. 

Once the image has been segmented, each granule needs to 
be identified. This operation, called connected-component label- 
ing, is an algorithmic application of the graph theory, where sub- 
sets of connected components are uniquely labeled based on a 
given heuristic. Connected-component labeling is used to detect 
the connected regions in the binary digital images produced by 
our segmentation algorithm. In this labeling process, a pixel be- 



longs to a granule if it shares at least one side with another pixel 
of the granule. 



3.2. Measuring the velocities 

Once the granules have been labeled in each frame, their tra- 
jectories are identified by comparing two consecutive images. 
During their lifetime, granules can split or merge into multi- 
ple objects. Disappearance or appearance between two succes- 
sive frames must be taken into account. Then, the life of coher- 
ent objects (i.e. granules) is defined between its appearance and 
disappearance if the granule does not split or merge. When the 
granule splits, the life of granule is stopped and its "children" 
are considered as new granules. In the same way, when gran- 
ules merge, the lives of the granules that merge are stopped and 
the new granule issued from the merging is considered as a new 
granule. Thus, we can follow a coherent structure between their 
birth and death. In the previous CST version, the temporal la- 
beling of the coherent structure was performed from one frame 
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Fig. 3. Full Sun Vx component from a two-hour sequence, spa- 
tial resolution 2.5 Mm. Arrows indicate the location of super- 
granule visible also on the Doppler map. 




Fig. 4. Full Sun Vy component from two-hour sequence, spatial 
resolution 2.5 Mm. Arrow indicates the location a supergranule 
visible also on the Doppler map. 



to the next by followi ng the bary center trajectory (iRoudier et al.l 
1 19991: iRieutord et al.l l2007l) and was essentially applied to the 
aligned data where large-scale motions were removed (solar ro- 
tation, drifts, etc.). However, when this method is applied to 
the HMI/SDO data, it generates some misidentifications of the 
barycenters between close granules because of the high value of 
the rotation velocity (around 2 km s"^). To avoid this extra noise, 
we now performed a time labeling of granules, that takes into ac- 



Fig. 5. Averaged Dopplergram from two-hour sequence where 
the solar rotation has been removed. Arrows indicate the location 
of supergranule. 
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Fig. 6. Solar grid where the dark circle indicates the limit of the 
validity of the horizontal velocity measurement. 



count a common granule area between image t and t -\- 1. This 
method ensures the temporal link between two snapshots of an 
evolving granule and improves their temporal labeling. 

Hence, each coherent structure (i) is defined by six values: 

1 . birth time Tbt and death time Tdi 

2. (Xi, Yi)h and (X/, Yi)d the positions at time Tbt and Tdi 

3. Vx and Vy are in the heliocentric-cartesian coordinates 
(lThompsonll2Q06h . respectively the horizontal component of 
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velocity in x {x increasing toward solar west) and in y (in- 
creasing toward solar north). 

From the birth and death locations and lifetimes of each co- 
herent structure, we derived a trajectory and mean velocity. To 
reduce the noise, coherent structures with a lifetime shorter than 
1 80 seconds were ignored. 

If one analyzes a long times series of images, it is useful to 
determine the time evolution of the velocity field; for this pur- 
pose a time window lasting 6t was used and trajectories were 
restricted to the time window. Hence, for a given time window, 
we derived a set of trajectories and velocities. The values of the 
velocities are of course not uniformly distributed in the field of 
view, and we need to know how they constrain the velocity field 
at a given scale: small-scale components are weakly constrained 
while large-scale ones are highly constrained. The maximum 
resolution for the velocity field is given by the density of trajec- 
tories. As pointed out in Rieut ord et al.l ('2001), granules cannot 
be used to trace plasma flow below a scale of 2.5 Mm (except 
for very rapid flows like in an "explosion" of granules); thus the 
determination of a large-scale flow needs a mesh size not smaller 
than 1.25 Mm. 

iRieutord et al.l (l2007l) indicated that the minimum size of 
the velocity mesh grid is related to the temporal resolution. 
They showed that one typical granule trajectory occupies a "vol- 
ume" of 1200Mm^s. When the "space-time" resolution does 
not reach this limit, many granules contribute to the veloc- 
ity in one mesh point. Their mean velocity is considered as 
the true local velocity but local fluctuations around this may 
also give some information on the local strength of convection. 
Because granules do not sample the field of view uniformly, 
a reconstruction of the velocity field along with its derivatives 
such as the divergence D = dxVx^ dyVy or the z-component 
of vorticity C = dxVy - dyVx requires some interpolation. Like 
iRieutord et al.l (|2007) we used the multi-resolution analysis for 
the interpolation, which limits the eff'ects of noise and error prop- 
agation. 

The smallest spatial and temporal resolutions achieved in the 
horizontal velocity map produced using the new CST algorithm 
are 2.5 Mm and 30 min, respectively. 

Fig. [2] shows some details around a sunspot in the intensity 
field and its segmentation where granules are visible; we plot an 
enlargement of the flow fields around the sunspot where diver- 
gent structures and the sunspot moat are clearly visible. 

On the full Sun white-light observation around 5x10^ gran- 
ules are detected. During a time tracking of 30 minutes around 
2.2 X 10^ coherent structures are identified. This number in- 
creases up to 7.5 X 10^ for a time sequence of two-hours. We 
show an example of the full Sun Vx and Vy in Fig. [3] and IH 
components from a two-hour sequence with a spatial resolution 
of 2.5 Mm. Arrows indicate the location of supergranules that 
are visible also on the Doppler map (Fig. O. Owing to the ra- 
dial flows within supergranules, two supergranules are marked 
in the Vx map close to the western limb where they are easier 
to identify than in the east- west component of the flow (and vice 
versa for the Vy map). The circle inside the solar grid in Fig. [6] 
represents the limit up to which the projected Vx and Vy are cor- 
rectly measured. Beyond this circle, granule tracking is difficult 
and errors increase rapidly as we near the limb. The useful field 
lies in the latitude and longitude range between ±60° in latitude 
and longitude for a two-hour sequence. Because this field covers 
a large part of the visible Sun area (80%), we define it as a quasi 
full-disk map of photo spheric horizontal velocities. 
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Fig. 7. Solar granulation from Hinode (1 pixel = 0^.^1089) on 
August, 30 2010 8h04mn33s (top left) and horizontal veloc- 
ities from a one-hour sequence (top right). Solar granulation 
from SDO (1 pixel = a.'5042 ), on August, 30 2010 8h04mn30s 
(middle left) and horizontal velocities from a one-hour sequence 
(middle right). Solar granulation from Hinode degraded to a 
pixel of a.' 5 on August, 30 2010 8h04mn30s (bottom left) and 
horizontal velocities from a one-hour sequence (bottom right). 



4, Precision 



Following iTkaczuk et al.l (l2007h . we can estimate the error prop- 
agation on the velocity of the CST algorithm. Over an image, 
granules are much less dense than pixels, therefore the veloc- 
ity field is sampled on a much coarser grid whose elements are 
of size 6 here 2.5 Mm. The velocity at a grid coordinate {x,y) 
is assumed to be the average of the velocity of granules whose 
average coordinates belong to the domain D around (x, 3;) (i.e. 
in [x - 612, X + ^/2], [3; - 6l2,y -h 612] ). The spatial resolution 
of the velocity field is given by the mesh size 6. If we assume 
that the dispersion on the barycenter is the same for the whole 
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Fig. 8. Velocity histograms from CST and LCT are plotted in 
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time series (cr), and that the time interval l^tk (of the granule k) is 
the same for all granules (i.e. all granules have the same lifetime 
AQ, then the ex pression of the dispersion of the average velocity 
(iTkaczuk et al.ll2007.) is 



dV = 



^fNAt' 



(1) 



where N is the number of trajectories falling in D. We found 
for the full Sun white-light SDO data for a time sequence of 30 
minutes and a mesh size of 2.5 Mm, that N=8. The precision of 
the location of the barycenter of granule is estimated to half a 
pixel, which represents 0V25 (183 km). With a median lifetime 
of 405 seconds SV = 0.22 km s"\ which is good enough at a 
scale of 2500 km, i.e. at the scale where granule m otion traces 
the large-scale plasma flows. iTkaczuk et al.l (l2007l) showed that 
precise velocity values need many granules in a grid element and 
a long-time interval. In other words, errors are less if a coarse 
resolution in space and time is used. 

When moving to the limb, the detection of granules becomes 
more difl&cult but the mesh size S covers a larger area on the 
Sun's surface. Both phenomena more or less cancel each other 
out and we find that close to the limb N=7 and SV = 0.24 km s"^ . 

Accordingly, we estimate that an error of SV = 0.25 km s"^ 
for the HMI/SDO data for a time sequence of 30 min is accept- 
able since most of the granules are showing velocities in the 
range 0.3-1 .8km s"^ in our analysis. This error seems to be large, 
but this is due in great part to the pixel size 0'.'5, which allows us 
to follow only the large granules on the Sun surface. However, 
when the velocities are averaged in space and time this error de- 
creases for example down to 0.06 km s"^ at the disk center with 
a window of 20 '' x 20 '' for a 24-h sequence (see below). This 
estimate error of the smoothed data is very similar to the one 
deduced from helioseismology analysis, which is found to be 30 
m/s for one day data (Table 2 of iSvanda et al.l (l201ll) ). 



5. Comparison of HINODE and SDO velocity fields 

To determine the quality of the velocities measured with white- 
light HMI/SDO data, quasi simultaneous observations of the 
Sun's surface with very high spatial resolution were obtained 
with the Hinode satellite during three-hours. The data were care- 
fully aligned, k-ci; filtered and resized at the same scale (1 pixel 
=0^.' 1089). The Hinode data were also degraded to the spatial res- 
olution of HMI/SDO for a detailed comparison. Three temporal 
sequences were treated exactly in the same way. We applied the 



CST to these sequences and show the results in Fig. [71 The cor- 
relation between velocities of the diff'erent sequences is about 
85% and one can easily recognize common patterns between the 
time series. The histograms of the velocity in Fig. [8] are very 
similar for all sequence, indicating that the amplitudes are also 
correctly measured. For comparison, local correlation tracking 
(LCT) was applied to the same one-hour sequence. This flow 
field is shown in the bottom of Fig. [71 Here we also find a good 
correlation (about 75%) between velocities computed by CST 
and LCT. The LCT velocities are lower, however than the CST 
velocities by 15%, likely becaus e of the spatial windo w size used 
in the LCT to convolve the data (iRoudier et al.ll 19991) . 

The high correlation between velocities and the quasi- 
identical amplitudes of the velocities indicate a good agreement 
between velocities measured with Hinode data (high spatial res- 
olution) and HMI/SDO data (low spatial resolution). This allows 
us to use HMI/SDO data with some confidence in determining 
the horizontal flows across the Sun. 



6. Power spectrum of the solar supergranulation 

The kinetic energy distribution among the scales is represented 
by the power spectrum of the velocity field. Like 'Rieut ord et al.l 
(2008), we computed the kinetic energy spectral density E{k) 
for a time window of 120 min and spatial field of 650 '' x 650 " 
centered on the disk center on August 30, 2010. 

The spectru m shown in F ig. [TOlis very similar to that exhib- 
ited in Fig. 2 of lRieutord"etal.l 12008). It exhibits the kinetic en- 
ergy contained in the supergranulation. With SDO data, we find 
the maximum of the spectral dens ity at a wavelength of 35.8 Mm 
quite close to 36 Mm found by iRieutord et al.l (l2008l) and the 
FWHM of the peak indicates that supergranulation occupies the 
rang e of scale s of [1 9, 56] Mm compared to [20, 75] Mm for 
■Rieutord et al.l (l2008l) . Our time sequence is shorter by a fac- 
tor 3.75 but our fiel d of view is 3.5 larger (400 '' x 300 '' for 
'Rieutord et al.ll2008l) . Like Rieutord et al., we indicate the best- 
fit power-law (3 and -2), which mimics the sides of the super- 
granulation spectral peak. At small scales below k = 0.lMm~^ 
we observe E(k) oc k due to the decorrelated random noise. Our 
present result confirms the previous findings. 



7. Solar differential and meridional flow 
determination 

There is a long history of inconsistent measurements of solar ro- 
tation (Beck 2000). One of the first applications of the CST algo- 
rithm on HMI/SDO data is to measure the horizontal velocities 
over 24h on the central meridian. Owing to computational time, 
we used a band of 50 '' along the central meridian and 48 time- 
sequences of 30 min. Then, 98 spatial windows of 50 '' x 20 '' 
were used to obtain the horizontal flow fields at diff'erent lati- 
tudes for each time- sequence of 30 min. Finally, we averaged all 
30 min time- sequences for each latitude to obtain the final Vx 
and Vy over 24h. The Vx are corrected for the Earth's orbital 
motion to obtain the sidereal rotation of the Sun, and Vy was 
also corrected for the BO evolution during 24h. This is the first 
time that solar rotation has been determined from granule dis- 
placement measurements over the full solar disk. Fig. |9] shows 
the diff'erential rotation computed in this way for four diff'erent 
dates close to the solar minimum. 

The standard deviation of the velocity field close to the disk 
center is 0.06 km s"^ (24h sequence). Owing to the change in 
size of the Sun in diameter of (1.4 pixel on the detector) (SDO 
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Fig. 9. Solar differential rotation on April 10, 2010, July 10, 2010, May 12, 2011, and May 13, 2011. We plot the rotation law 
determined by spectroscopic method (' Howard & Harveylll970i) in the four panels oj 
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Fig. 10. Kinetic energy spectra of the horizontal velocities. The 
vertical dashed line emphasizes the 10 Mm scale, which is usu- 
ally taken as upper limit of mesogranular scale. Two power laws 
are shown on each side of the peak, as well as that of the small- 
scale noise. 



Earth orbit) and also to the BO evolution during 24h, the field 
of good measurement was reduced to -50° and -\- 45° in lon- 
gitude and ±50° in latitude around the disk center. This limita- 



tion will be improved in the near future (new observations on 10 
December 2011) by reducing the eff'ects of dilatation and con- 
traction of the Sun's diameter on the CCD. 

Our results of April 10, 2010 agree with the spectroscopic 
determination of the solar rotation of Howard & Harvey (1 970h 
but the rotation seems slightly faster on the other dates with an 
equatorial velocity of 1.99 km s"^. The daily Wolf numbers are 
8, 14, 33, an d 26 and 8, 16.1, and 41.6 monthly. We expect dif- 
ferences with lHoward & Harvey! (Il970h . simply because we are 
not studying the same time period. The solar diff'erential rotation 
appears to change in time particularly at high latitudes, but we 
need a more extensive analysis to confirm these variations. 

The meridional flows with Vrms = 90 m s"^ given by the Vy 
component for 24h, is plotted in Fig.[TT]for a quiet Sun on April 
10 2010 between -50° and -h 45° in latitude. Because the Vrms 
of the mean meridional flow is 90 m s"^ the expected signal of 
20 m s"^ is completely hidden. Several days are required to re- 
duce the Vrms amplitude and determine the meridional velocity. 
This work is scheduled for the near future because we must also 
include the Sun's diameter evolution on the SDO detector over a 
long period of time (Earth orbit). 

8. Discussion and conclusions 

Determining horizontal velocity fields on the solar surface is cru- 
cial for understanding the dynamics of the phostosphere, the dis- 
tribution of magnetic fields and its influence on the structures of 
the solar atmosphere (filaments, jets, etc.). 



Roudier et al.: Quasi full-disk maps of solar horizontal velocities using SDO/HMI data. 
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Fig. 11. Meridionnal velocity at different latitudes across the 
central meridian on April 10, 2010, quiet Sun. 



We have shown that the reconstruction of the velocity field 
can be achieved over the Sun's entire visible surface with 
HMI/data by using CST. From this reconstruction, by using the 
multi-resolution analysis, one can also compute the velocity field 
at diff'erent scales along with its derivatives such as divergence 
and curl. An uncertainty in the velocity of about 0.25 km s"^ 
for a time sequence of 30 min and a mesh of 2.5 Mm is accept- 
able if compared to the granule velocities which range between 
0.3 km s"^ and 1.8 km s"^ The kinetic energy power spectra 
obtained on a large field of the Sun confirms the previous quan- 
titative studies of the dynamic of the Sun's surface. The first de- 
termination of the solar rotation using granule tracking demon- 
strates the potential of our method for studying the evolution of 
solar rotation on short- time scales. The method can also be used 
to remove the solar rotation for studying smaller scale evolution 
(supergranulation flows, etc.). 

In a near future, it will be possible to derive the velocities 
on the solar surface in the local coordinates system (V^, Vg, V^) 
from the line-of- sight component velocity (Vdop) and horizon- 
tal velocity fields on the solar surface (Vx,Vy). In particular, it 
will be possible to compare the spherical analysis of the Sun's 
motions to those of numerical simulations, thus improving our 
knowledge of the physical properties of the Sun. 

Contemplating the CST method and the HMI/SDO data, 
many applications can be envisioned on time scales ranging 
from 30 minutes to several days or even months with spatial 
scales ranging form 2.5 Mm up to the radius of the Sun. This 
opens a new field of study allowing, e.g., 3. measure of the influ- 
ence of photospheric motions on the dynamics of the outer so- 
lar atmosphere such as the triggering of filament eruptions, see 
iRoudier et al.l (l2008h . Many other applications will be possible 
on phenomena that surface motions involve. 

Today, 12 min of processor time are required to compute 
30 min of solar observations. Improvement of the code is in 
progress. Currently, the Fortran90 CST code is parallelized in 
a shared memory approach involving OpenMP. We experienced 
good scaling on a production data set (full sun) with 73.16% 
parallel efficiency on 24 cores. Tests were conducted on an SGI 
Altix UV ccNUMA node with 96 cores and 1TB of RAM. This 
CST code (Fortran90) will be implemented in the German Data 



Center for SDO and in the Joint Science Operations Center, and 
will be available to solar community. 
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